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We introduce a multiplex particle focusing phenomenon that arises from the hydrodynamic interaction 
between the viscoelastic force and the Dean drag force in a microfluidic device. In a confined microchannel, 
the first normal stress difference of viscoelastic fluids results in a lateral migration of suspended particles. 
Such a viscoelastic force was harnessed to focus different sized particles in the middle of a microchannel, and 
spiral channel geometry was also considered in order to take advantage of the counteracting force, Dean drag 
force that induces particle migration in the outward direction. For theoretical understanding, we performed 
a numerical analysis of viscoelastic fluids in the spiral microfluidic channel. From these results, a concept of 
the £ Dean-coupled Elasto-inertial Focusing band (DEF)' was proposed. This study provides in-depth 
physical insight into the multiplex focusing of particles that can open a new venue for microfluidic particle 
dynamics for a concrete high throughput platform at microscale. 

Viscoelastic fluids, which can be encountered extensively in nature, demonstrate extraordinary physical 
phenomena due to their intrinsically viscous and elastic features. There are some intriguing examples 
related to the viscoelasticity: most biological cells produce extracellular polymeric substances (EPS), such 
as polysaccharides, nucleic acids, polypeptides, proteins, growth regulators, lipids, and a number of other com- 
pounds not yet fully characterized 1 " 4 . The protein secretion generates natural viscoelastic substances to protect 
cells from the harmful effect of toxic substances 3 ' 5 . The most common fluid in our daily lives is saliva, the 
viscoelasticity of which forms a distinct beads-on-a-string (BOAS) structure 6 . Lastly, cells in a Newtonian 
medium experience viscoelastic deformation during the flow throughout vessels provoking viscoelasticity- 
induced force, which is used to isolate rare target cells from heterogeneous sample 7,8 and to separate normal 
erythrocytes from malaria-infected erythrocytes 910 . 

Very recently it has been revealed that the positive first normal stress difference (iVj) generated by the 
viscoelastic nature leads suspended particles or cells to migrate into certain positions in a confined microchan- 
nel 11 " 18 . Manipulating particle dynamics with help of the intrinsic properties of the medium is very desirable for 
easy, cost-effective, and continuous separation compared to conventional techniques such as filtration 19 " 21 , 
electrophoresis 22 , dielectrophoresis 23 , and hydrodynamic chromatography 24 " 26 . Worth noting is that particles 
in the viscoelastic Poiseuille flow are migrated towards the center and four corners of the channel, i.e., the position 
with the minimum shear rate 1318 . Moreover, the number of multiple equilibrium positions can be reduced to one 
equilibrium at the centerline by increasing the flow rate due to the synergetic effect of inertia and viscoelasticity 18 . 
The lateral migration of particles in the viscoelastic fluid is driven by the imbalance of the compressive nonlinear 
elastic forces, which are described in terms of the first and second normal stress differences, Ni (y) = o xx — a yy and 

N 2 (y) = (Tyy-(Tzz. 

The Dean vortices, one of the secondary forces, can also contribute to particle migration. The fluid within the 
curved channel should satisfy the mass balance between inner-wall and outer-wall regions. This induces a 
pressure gradient in a radial direction, thereby generating two counter rotating vortices at the channel. This 
secondary flow is characterized by the Dean number, De = Re c -\/ Dh / 2R, where Re c is the channel Reynolds 
number, D h is the hydrodynamic channel height, and R is the radius of curvature of the spiral channel. A recent 
study presents the balance between the inertial lift force, F L , and the Dean drag force, F D ~ pU^aD^r -1 , which 
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Figure 1 | Schematic representation of the Dean-coupled Elasto-inertial Focusing. Different sized particles are randomly injected into an inlet and 
displaced with differential lateral migration at an outlet. In a confined spiral channel, particles experience both elastic force and Dean drag force, the 
interaction of which determines equilibrium positions of the suspended particles. 



determines the preferred location of particles in a curved channel 27 ' 28 . 
When the density of a particle is larger than that of the surrounding 
medium, the particle subsides. In a similar way, larger particles can be 
collected in an outer wall with the aid of the characteristic velocity 
distribution on the cross-sectional plane normal to the main flow 29 . 
As the density of a particle is similar to that of the medium, the ratio 
of F L /F D is a determinant of the equilibrium positions. For instance, 
larger particles with higher F L /F D have a lateral migration towards 
the inner wall 30 ' 31 . Thus, the scaling law between F L and F D acts as a 
dominant focusing factor: if the Dean flow is dominant (i.e., F D 
> F L ), there is no focusing. Otherwise, focusing is favored due to 
the dominant inertial lift force. The scaling law, F L /F D ~(l/d) 
(a/Dh) 3 Re" («<0) containing parameters of Re c and curvature 
ratio, S, indicates that the ratio of particle size and channel dimen- 
sion, alD h , plays an important role in determining the focusing 
behavior. For instance, small particles with large channel dimensions 
would not be focused. Preferential focusing of particles can only be 
achieved when a/D h is larger than 0.07 30 ' 32 ' 33 . The aforementioned 
inertial microfluidics demonstrates that the defocusing of a particle 
occurs as F D and Re c increase. However, an equilibrium state is 
unstable since the superposition of inertial lift and Dean drag forces 
gives two modified stable equilibrium positions at which the balance 
is susceptible to be broken by the Dean vortices 28 . Consequently, no 
equilibrium position exists in that a particle size is much smaller than 
a microchannel or the balance of F L /F D breaks the criterion. In this 
respect geometry itself does not necessarily guarantee complete par- 
ticle manipulation in a microchannel, while the elastic force driven 
by fluidic nature (i.e., the differential of Nj) can solely focus sus- 
pended particles in a line stably 715 . 

Here, we explore an extraordinary particle focusing in a micro- 
fluidic channel by probing two counteracting forces: one is an elastic 
force induced by the nature of a viscoelastic medium and the other is 
a Dean drag force resulting from the curved channel geometry. Also, 
we investigate the physical reasons for the multiplex particle focusing 



by numerical analysis and propose a new concept of the 'Dean- 
coupled Elasto-inertial Focusing (DEF)' band. 

Results 

Different sized particles are randomly injected into an inlet of a spiral 
channel using a viscoelastic medium (Fig. 1). In the confined spiral 
channel the particles undergo both of elastic force (F E ) and Dean 
drag force (F D ), which determines equilibrium positions of the sus- 
pended particles. As a result, The DEF band is created on the right 
side of the confined spiral channel (see the red-box in the cross- 
section of the channel). A narrow spiral channel with an aspect ratio 
(y) is 0.25 shows different distribution of the first normal stress 
difference induced by the elastic nature of viscoelastic fluid, resulting 
in the focusing band and displacement of the different sized particles 
as shown in the outlet of the schematic representation. It should be 
noted that once particles are focused in the region of the DEF band, 
lateral migration is determined by the competition of F D and F E . 

No particle focusing is observed at the outlet when each polysty- 
rene (PS) particle with sizes of 1.5 um, 5 um, and 10 urn is used in a 
Newtonian medium (Fig. 2a ~ c and SI Fig. 1). However, particles 
were successfully focused when using a viscoelastic medium. Single- 
line particle focusing has been successfully achieved by using a vis- 
coelastic medium, which is attributed to the combination of F D and 
F E . Also, the larger particles were found to move closer to the outer 
wall than the smaller particles. We also examined multiplex particle 
focusing using both small and large particles with diameters of 
1.5 um and 10 um in a viscoelastic medium (Fig. 2g). It is expected 
that the particle to particle interaction would increase with increasing 
particle concentration, which induces the perturbation of flow and 
then prohibits the particle focusing. In this sense, a microfluidic 
system using hydrodynamics has not been regarded as a right tool 
for focusing size -dependent multiple cells or particles in a robust 
manner. It is shown that both the small and large particles flow 
randomly at the inlet and the randomness is still present at the 
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Figure 2 | Experimental results of particle dynamics in both the Newtonian and viscoelastic mediums, (a) , (b) and (c) PS particles are not focused in the 
Newtonian medium but focused in the viscoelastic medium, (d) Dual PS particles with 1.5 and 10 (am are not focused in a Newtonian medium, 
(e) Probability distribution function of 1 .5 urn PS particles with respect to the number of turns, (f) Probability distribution function of 10 (am PS particles 
with respect to the number of turns, (g) Dual PS particles with 1.5 and 10 (am are successfully focused with differential lateral migration in the viscoelastic 
medium, (h) Probability distribution function of 1.5 (am PS particles with respect to the number of turns, showing particle focusing at the center, 
(i) Probability distribution function of 10 (am PS particles with respect to the number of turns, showing particle focusing near the outer wall. 



second turn of the spiral channel (Fig. 2g). After the particles experi- 
ence the tenth turn, both the particles have different equilibrium 
positions, indicating that 1.5 (am PS particles are at the center line 
and 10 (am PS ones are close to the outer wall (Fig. 2g). The prob- 
ability distribution function of both 1.5 and 10 (am PS with differ- 
ential lateral migration elucidates the preferential particle focusing at 
the center for the small particle and near the outer wall for the large 
particle (Fig. 2h and 2i). However, both particles are not focused in a 
Newtonian medium even though the particles go through the tenth 
turn (Fig. 2d). The probability distribution function of both 1.5 and 
10 (am PS with differential lateral migration shows that no particle 
focusing occurs regardless of the particle size (Fig. 2e and 2f). The 
results demonstrate that multiplex particle focusing is achievable 
using a hydrodynamic system through a viscoelastic medium with- 
out any external forces. 

The different sized particles are entrained parallel, and bifurcate at 
the outlet: 1.5 (am PS particles towards the upper outlet and 10 (am 
PS ones towards the lower outlet (Fig. 3 and SI movie 4). The large 
particles with the diameter of 10 (am move closer to the outer wall of 
the microchannel due to the combination of the elastic and Dean 
drag forces. However, the small particles with the diameter of 1.5 (am 



are positioned at the center because the effect of the Dean drag force 
acting on the particles is comparably smaller than that of the elastic 
force. The lateral displacement of both particles was acquired and 
fitted by the Gaussian distribution function (see the inset in Fig. 3). 
There is no overlap between the lateral displacements of 1.5 (am and 
10 (am PS particles, which guarantees high separating efficiency at 
the outlet. 

It turns out that the high magnitude of AT? appears near the inner 
wall in the cross-section of a straight channel (Fig. 4a and 4c). Nj is 
not found in a Newtonian medium (meaning Wi = 0) since there is 
no elastic nature in the Newtonian fluid, thus being subject to only 
the Dean drag force and inertial lift force. The minimum values are 
found at four corners and the center as Wi becomes large, which 
leads particles to the center and the phenomenon is called as the 
'elasto-inertial effect' 18 . The value of iVl/liV]™ 1 *! shows symmetry 
along the lateral position of the channel and increases with increasing 
Wi (Fig. 4c). In addition, the curvature of the spiral channel creates 
Dean vortices causing the lateral migration of the flow in the channel. 
Therefore, larger particles undergo more lateral migration on 
account of the increased De p than smaller particles (Fig. 4b and SI 
movie 2). De p corresponds to the particle dean number with 1.5, 5, 
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and 10 um PS, respectively. It turns out that the high magnitude of 
iVj appears near the inner wall in a spiral channel not like the four- 
fold symmetry in the straight channel (Fig. 4d and SI movie 3). It 
shows that VN\ j | VAT™ | is an indicator of the elastic force exerted 
on the particle because of F E ~ a 3 VNi, which represents particles can 
migrate toward the outer wall (Fig. 4e). This asymmetric gradient of 
iVj leads the equilibrium position to be off-centered. 

Discussion 

An elastic force, F E , can be introduced when a viscoelastic fluid is 
used as a medium 13,34 " 36 . The elasticity, inherently present in dilute 
polymer solutions, can yield particle focusing into the mid-plane of 
the channel. The particles have different spatial distribution depend- 
ing on their sizes and rheological properties of the viscoelastic med- 
ium. The primary nonlinear elastic force affecting the lateral 
migration is N ly since the value of N 2 is normally much greater than 
N 2 15,18,37 . The relative magnitude of elastic and inertial forces is 
assessed by the Elasticity number, El = Wi/Re 38 . Nj is associated 
with Wi = ky 0 which is zero when the fluid is purely viscous, i.e., no 
elasticity, and increases with increasing the elasticity of fluid. It is 
worthy to note that Wi is represented as a physical factor to dem- 
onstrate the elasticity of the viscoelastic fluid since Wi is determined 
by the characteristic relaxation time of polymer, a physical quantity 
of time required for a perturbed polymer chain to return to the 
equilibrium. Therefore, it is important to characterize Wi to find 
out an optimal condition for particle focusing from the Wi-De dia- 
gram. The use of the elastic force is indeed an appealing method to 



focus particles at the center of the channel without any external 
forces and complex geometry. In the literature, small particles with 
(a/D h ) < 0.04 may not successfully be focused due to insufficient 
length of a straight microchannel 18 . The present spiral channel has 
approximately 50 cm of arc-length and the length is 10 to 25 times 
longer than the straight microchannel, which creates the effective 
elastic force and Dean drag force for the sufficient particle migration 
regardless of the particle size, e.g., even under (a/D h ) < 0.04. 

The elastic force generates equilibrium positions at the center and 
four corners due to AT? acting on the particle 818 ' 38 (Fig. 4a and 4c, and 
SI movie 1). Particles can be focused at the center as inertial force is 
increased by increasing the flow rate and F E becomes large enough 
(meaning, Re > 0 and Wi > 0). On the other hand, F D increases with 
an increase in Re c according to the following equations: De = 
Re c -^/Dh/2R and F D ~ func(a,De) 31 . Since F D is a function of the 
particle size and De, the lateral migration behaviors of particles are 
distinct (Fig. 4b). Here, another Reynolds number is adopted 28 as the 
particle Reynolds number, Re p = Re c (a 2 /Dfy = U m a 2 /vDh, includ- 
ing parameters for the perturbed channel flow due to suspended 
particles. We claim that the suspended particles are also perturbed 
by the secondary flow in a curved spiral channel so that the particle 
De can also be given by De p = Re p ^Jd^/IR = Re c (aj Dh ) 2 (Dh / 2R) . 
Hence, the Dean drag force changes significantly with respect to a 
particle diameter, thereby resulting in a differential lateral migration 
which is essential for the differential particle focusing. As a con- 
sequence, particles have equilibrium positions where the elastic force 
is determined by De P . 
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Figure 4 | Numerical simulation of the Dean flow and the first normal stress difference, (a) Results of the first normal stress difference with 
Wi = 0, 0.4, 0.7, and 1 in a straight channel with a channel aspect ratio of 1. The minimum values are found at four corners and the center, (b) Results of 
the Dean flows for De = 0.22, 4.33, 13, and 21.6 in a curved spiral channel with a channel aspect ratio of 1. De p corresponds to the particle dean number 
with 1.5, 5, and 10 (am PS, respectively, (c) Normalized first normal stress difference in a straight channel with a channel aspect ratio of 1. It shows 
nearly zero at Wi = 0 and increases with increasing Wi. (d) Results of the first normal stress difference for Wi = 0, 0.5, and 1 in a curved spiral channel with 
a channel aspect ratio of 0.25. It shows an asymmetry, resulting in a lateral migration of particles towards the outer wall, (e) Normalized differential value 
of the first normal stress difference in a curved channel with a channel aspect ratio of 0.25. The differential value of the first normal stress difference is an 
indicator of the elastic force exerted on the particle as given by the equation, F E ~ a 3 VNi. 



The main forces related to the present system acting on the sus- 
pended particles are F D and F E . Here, De P and F D ~6ninavDe P /Dh 
are equated to derive F D = 6rjvRe 2 c S(a/Dh) 5 , and F E is expressed as F E 
= Aa 3 WNi, which is proportional to the particle size and the differ- 
ential of Ni. The two forces lead to a scaling law of (F D /F E )~ 
a 2 3rjRe 2 c /AVNiD 5 h , which helps one to understand physical para- 
meters associated with the equilibrium positions of particles with 
distinct sizes. F D and F E acting on the particles in the microfluidic 
system were evaluated with use of the calculated quantities, VA/i, 
U De an> R e c> an d De P . The relevant results are presented in Fig. 5. 
F D significantly increases with decreasing channel height, which 
implies that the suspended particles tend to move in an outer dir- 
ection of the channel. The normalized Dean forces significantly 
increase as the channel aspect ratio (y) decreases, thereby inducing 
particle-size dependent lateral migration (Fig. 5a). That is, the par- 
ticle size is one of the key factors to govern the multiplex particle 
focusing phenomenon. Conflicting results may be found in two 
papers: first, Carlo et al. 27 demonstrated that F D decreases with 
increasing particle size in a given channel aspect ratio since the lateral 
force and Dean force are complementarily counteracted as 
(F L /F D )~(l/d)(a/D h ) 3 R*. Second, Papautsky et al. 31 addressed 



the dependence of F L /F D on the particle size, namely the decrease 
of F D with respect to the particle size. Hence, as the ratio increases, 
the lift force mainly acts on the particles. Since the two researches, 
however, adopted Newtonian mediums in their experiments, the 
mid-plane particle focusing band caused by the elastic nature of a 
viscoelastic medium could not be formed. Therefore, the Dean flow 
effect on the particle migration could not be developed. In this 
respect, the proposed multiplex particle focusing offers a noteworthy 
advantage on differential migration of different size particles as a 
result of the competition between F E and F D . The former is created 
by the elasticity of the viscoelastic fluid, and the latter is generated in a 
focusing band region by the geometry of the curved spiral channel. 

In addition to the aforementioned, decreasing y is advantageous 
for the reduction of Focusing Band Ratio (FBR) which is determined 
by the velocity distribution of Dean flow (SI Fig. 2). Here, FBR is 
defined as (Focusing Band Height/ Channel Height) X 100 [%], which 
means that the particle-focusing-band-height also decreases with a 
decrease in FBR. The schematic illustration is presented in SI Fig. 2b. 
The focusing band heights are determined by outward velocity dir- 
ection in the mid-plane region where particles are focused in the 
mid-plane (SI Fig. 2a). FBR reaches 44% when y is below 0.5, indi- 
cating that the band height is comparable to a particle size of 1.5 |im 
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Figure 5 | Numerical results of the Dean drag and elastic forces at the horizontal centerline of the microchannel when Wi = 1. (a) Normalized Dean 
force with respect to channel aspect ratio (H/W) and particle size ranging from 1.5 to 10 (am. A boxed area indicates possible conditions for inducing 
particle migration, (b) Dean drag forces are normalized with the maximum Dean drag force, indicating the effect of the ratio of particle size to the 
hydraulic diameter of a microchannel. An inset graph shows both of the dean drag and elastic forces normalized by the maximum Dean drag force under 
the channel aspect ratio of 0.25. (c) Normalized lateral equilibrium positions of PS particles with 1.5, 5, and 10 (am diameters. A schematic illustration 
represents net forces acting on the different-sized PS particles. 



which results in successful focusing in the mid-plane particle focus- 
ing and lateral migration (SI Fig. 2b). The ratio of the particle size to 
the microchannel dimension, a/D h , determines whether the system 
leads to a differential particle migration due to the interaction 
between F D and F E or not. The Dean forces remarkably increase with 
respect to y when a/D h is larger than 0.06 (Fig. 5b), and are almost the 
same when y is 0.5 and 0.25. These results reveal that y = 0.25 is 
sufficient to provoke the differential lateral migration by using F E 
for the mid-plane focusing and F D for the lateral migration. 
Furthermore, the magnitudes of F D and F E are significantly different 
as alD h becomes larger than 0.08 in an inset graph; e.g., our system 
with 100 |im in width and 25 |im in height shows that particle 
migration is not likely to occur since the magnitudes of F D and F E 
are almost the same causing no difference in net forces acting on the 
particle. As the particle size is larger than 5 (im, meaning a/D h is 
greater than 0.08, F E becomes more significant than F D , which leads 
to the differential particle migration. 

The equilibrium positions of particles are determined at which the 
magnitudes of F D and F E are identical. Despite the fact that it is 
difficult to obtain the exact magnitudes of F D and F E quantitatively, 
the equilibrium positions according to the different particle sizes can 
be analyzed by using the scaling law. When the particle size is around 
1.5 (im, F E is almost the same as F D . Therefore, the particles are not 



laterally migrated, but are focused at the center (Fig. 5c). When the 
particle size is large enough, e.g., the size is as large as 10 (im, F D 
becomes greater than F E . Thus, the equilibrium position is located at 
a distance far from the center. This extraordinary particle dynamics 
according to the particle size is explained with the combination of the 
elastic force inducing the mid-plane particle focusing and the Dean 
drag force leading to the particle defocusing from the mid-plane. 
That is, the combinational manipulation of F D and F E is essential 
to outperform the ability of particle focusing and separation. F D and 
F E play an important role for the separation of different sized part- 
icles simultaneously. 

The viscoelastic moduli G' and G" were measured with respect to 
the angular frequency, which gives the characteristic relaxation time 
for different kinds of PEO solutions assuming the Maxwell model 
described by G'(co) = f]lco 2 /(l + l 2 w 2 ) andG"(co) = rjw/(l + Pco 2 ), 
where r\ is the viscosity, co is the angular frequency, and X is the 
characteristic relaxation time. The characteristic relaxation times 
for 500, 1000, 3000, and 5000 ppm PEO solutions are estimated to 
be 0.539 X 10 3 , 1.559 X 10 3 , 3.522 X 10 3 , and 16.1 X 10 3 , respectively. 
The flow rate was increased from 50 to 750 (il/h to evaluate the value 
of Wi = XQ/hw 2 and plot the Wi-De diagram, where h is the height 
and w is the width of the curved spiral channel. It is found that 
multiplex particle focusing initiates when Wi becomes larger than 
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Figure 6 | Experimental results for particle focusing as a function of De 

and Wi. No focusing corresponds to red-crosses, single-line focusing blue- 
triangles, and multiplex focusing pink-circles. An inset graph indicates a 
separation efficiency with respect to Wi and De. 

2 and De is greater than 10 (Fig. 6). The separation efficiency of the 
proposed system was plotted with respect to De (the inset in Fig. 6). 
The efficiency is very low when De ~ 2.4. This is because particles 
only undergo the elastic force within the channel and are not effec- 
tively focused at equilibrium positions depending on their sizes. As 
the Dean drag force becomes large as much as De ~ 25, the efficiency 
becomes extremely high compared with De ~ 2.4 and De ~ 18, the 
efficiency of which is comparable to the highest separation efficiency 
reported in the literature 31 ' 38 " 42 . This microfluidic system is capable of 
acting as a versatile platform that requires only the viscoelastic char- 
acteristics of a medium without any additional external forces. 

In conclusion, we assert that multiplex particle focusing can be 
achieved by the fine tuning of the two main forces: the elastic force 
driven by the nature of a viscoelastic medium and the Dean drag 
force induced by a curved spiral channel. The particles are situated at 
an equilibrium position between the two main forces. A numerical 
simulation was carried out to fully understand the multiplex particle 
focusing phenomenon. The particle migration towards the outer wall 
in a focusing band is increased as the particle size becomes larger. 
Also, it is found that the channel aspect ratio is important to focus 
particles into the mid-plane of the channel, and that when the focus- 
ing band dimension can reach the particle size the particle movement 
is not perturbed by the Dean vortices but influenced by the lateral 
Dean drag force and elastic force. It is noticed that particle size- 
driven multiplex focusing is possible through hydrodynamics 
encompassing the elastic force and the Dean force. We expect that 
this study will be able to provide a comprehensive understanding on 
the particle focusing mechanism which is applicable to various 
industrial applications ranging from labeling, particle and cell sepa- 
ration, and isolating functionalized particles for biomedical and 
pharmaceutical practice. 

Methods 

Experimental methods. Three types of PS particles with 1.5, 5, and 10 um diameters 
were used in this study. 1.5, 5, and 10 um PS particles were purchased from the Bangs 
Laboratories. A Newtonian fluid was prepared by mixing 22 wt% Glycerin with 
distilled water to match the densities between the medium and the PS particles. For a 
viscoelastic fluid, aqueous solutions of polyethylene oxide (PEO, M w = 1,000,000) 
with 500, 1000, 3000 and 5000 ppm were used. The Flow rate was increased from 50 
to 750 |il/h so as to evaluate various Wi and De. The PS particles were suspended in 
both the Newtonian and viscoelastic solutions with a small amount of surfactant 
(Tween 20, Sigma- Aldrich). The solutions were sonicated for 5 minutes for better 
dispersion of each particle. Spiral channels with the initial radii of curvatures of 4, 6.6, 
8.25, and 9.5 mm were fabricated to evaluate different De. All the spiral channels have 



ten loops, total arc-lengths are approximately 50 cm, and the width and height are set 
to 100 and 25 um. 

Numerical simulation. We carried out a finite element simulation using the 
Oldroyd-B model in straight and curved microchannels to solve the viscoelastic 
behavior of the fluid. For the simulation, governing equations are non- 
dimensionalized by using Re c and Wi. For steady state, the momentum equation is 
expressed as Re(u-V)u = V-( — pi + (n s /r])[Vu + (Vw) r ] + T) and the extra stress 

contribution by the viscoelasticity becomes T+WiT = (jl p /^ + (Vw) r j , 

v v dT 

where T is the upper convected derivative operator expressed as T = + 

(wV)T- J(Vw)''T + r-(Vw) r j , X, is the characteristic relaxation time, rj s is the relative 

polymer viscosity, n s is the relative solvent viscosity, and the total viscosity is n = r\ s + 
rj p . The two equations along with the continuity equation are coupled to solve the 
viscoelastic fluid by varying Wi in order to observe the effect of Nj in the straight and 
curved microchannels. The relative solvent and polymer viscosities were set to 0.6 
and 0.4. 
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